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We derive the perturbation equations for relativistic stars in scalar-tensor theories of gravity 
and study the corresponding oscillation spectrum. We show that the frequency of the emitted 
gravitational waves is shifted proportionally to the scalar field strength. Scalar waves which might 
be produced from such oscillations can be a unique probe for the theory, but their detectability is 
questionable if the radiated energy is small. However we show that there is no need for a direct 
1^-^ ■ observation of scalar waves: the shift in the gravitational wave spectrum could unambiguously signal 

' the presence of a scalar field. 
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Cn 1 I. INTRODUCTION 

<n : 

| Scalar-tensor theories of gravity are an alternative or generalization of Einstein's theory of gravity, where in addition 
. to the tensor field a scalar field is present. The theory has been proposed in its earlier form about half century ago 
ylclla, and it is a viable theory of gravity for a specific range of the coupling strength of the scalar field to gravity 
■ QLuLul- Actually, the existence of scalar fields is crucial (e.g. in inflationary and quintessence scenarios) to explain 
the accelerated expansion phases of the universe. In addition, scalar-tensor theories of gravity can be obtained from 
the low-energy limit of string theory and/or other gauge theories. Experimentally, the existence of a scalar field has 
not yet been probed, but a number of experiments in the weak field limit of general relativity set severe limits on the 
existence and coupling strengths of scalar fields [|| Q • 
— «. . A basic assumption is that the scalar and gravitational fields ip and g^ are coupled to matter via an "effective 
metric" <7 M „ = A 2 ((p)g^ l/ . The Fierz-Jordan-Brans-Dicke QH^EI theory assumes that the "coupling function" has the 
I , form A(ip) — aotp, i.e., it is characterized by a unique free parameter — (2^bd +3) _1 , and all its predictions differ 
h~L ' from those of general relativity by quantities of order «g [g . Solar system experiments set strict limits in the value 
. . | of the Brans-Dicke parameter wbd, i.e., wbd ^ 40000, which suggests a very small cvq < 10~ 5 (see 0,0)- 

In the early 1990s, based on a simplified version of scalar tensor theory where A(tp) = a$ip + f3tp 2 /2, Damour and 
Esposito-Farese H,E3 found that for certain values of the coupling parameter (3 the stellar models develop some strong 
field effects which induce significant deviations from general relativity. This sudden deviation from general relativity 
for specific values of the coupling constants has been named "spontaneous scalarization" . Harada studied in 
more detail models of non-rotating neutron stars in the framework of the scalar-tensor theory and he reported that 
"spontaneous scalarization" is possible for (3 < —4.35. DeDeo and Psaltis suggested that the effects of scalar fields 
might be apparent in the observed redshifted lines of X-rays and 7-rays observed by Chandra and XMM-Newton [T^ ] 
and in quasi-periodic oscillations (QPOs) |l3j. 

Recently, we have examined the possibility to obtain the information for the presence of the scalar field via grav- 
itational wave observations of oscillating neutron stars (0], hereafter Paper I). This previous study has been done 
using the so-called Cowling approximation. In this approximation one studies the fluid oscillations freezing the per- 
turbations of the spacetime and of the scalar field. Even under these assumptions the effect of the scalar field on 
the fluid perturbations can be significant. We showed that for values of < —4.35 the oscillation frequencies of the 
fluid change drastically, and the observation of such oscillations can verify or rule out the spontaneous scalarization 
phenomenon. 

It has been suggested that stellar oscillations can provide a unique tool for estimating the parameters of the star, 
i.e., mass, radius, rotation rate, magnetic field and equation of state. These ideas have been developed in the last 
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decade in a series of papers [H E IH El El IH E|, where the properties of the various families of oscillation 
modes have been used to probe the stellar parameters. The modes which are mainly excited during the formation of 
a neutron star or during starquakes and emit detectable gravitational waves are the fluid / and p modes and the w 
modes, which are associated to oscillations of the spacetime |22j| . 

The effect of the scalar field on the / and p modes has been examined in Paper I (in the Cowling approximation). 
In this article we derive the full set of equations describing the oscillations of a relativistic star, i.e., the perturbations 
of the fluid, the spacetime, and the scalar field. Since the stellar models are spherically symmetric, the oscillations 
can be classified as axial or polar depending on their parity, and we can derive perturbation equations for each class 
of perturbations. In the polar case we show that the wave equations describing the perturbations of the fluid and 
spacetime couple to the wave equation describing the perturbations of the scalar field. In other words the polar 
perturbations are affected not only by the presence of the scalar field in the background but also by the coupling 
with the wave equation describing the perturbations of the scalar field. In the axial case we find a single equation 
describing the perturbation of the spacetime. This equation is not coupled to perturbations of the fluid and of the 
scalar field: the scalar field only affects the dynamics through its influence on the background. 

The paper is organized as follows. In the next section we establish our notation and briefly introduce the theoretical 
framework for the scalar-tensor theories of gravity we consider. In Section IlIII we derive the perturbation equations 
which will be used for the numerical estimation of the oscillation frequencies. In Section Hvl we describe the methods 
that have been used to derive the axial w modes in scalar-tensor theory of gravity and show the results. In the final 
Section we discuss the results and their implications. We have included in three Appendices the details of the 
various analytic and numerical calculations. In Appendix 1X1 we describe the perturbations of the energy momentum 
tensors for the fluid and the scalar field, while in the next Appendix[B]we provide the analytic forms of the perturbed 
Einstein equations. Finally, in the last Appendix [U] we describe the numerical techniques that have been used to 
calculated the quasinormal modes. In this paper we adopt the unit of c = G = 1 , where c and G denote the speed of 
light and the gravitational constant, respectively, and the metric signature of (— , +, +, +). 



II. STELLAR MODELS IN SCALAR-TENSOR THEORIES OF GRAVITY 

In this section we will study neutron star models in scalar-tensor theory of gravity with one scalar field. This is 
a natural extensions of Einstein's theory, in which gravity is mediated not only by a second rank tensor (the metric 
tensor g^) but also by a massless long-range scalar field <p. The action is given by Q 
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where all quantities with asterisks are related to the "Einstein metric" g^ u , then i?* is the curvature scalar for 
this metric and G* is the bare gravitational coupling constant. \& TO represents collectively all matter fields, and S m 
denotes the action of the matter represented by \& m , which is coupled to the "Jordan-Fierz metric tensor" g^. The 
field equations are usually written in terms of the "Einstein metric" , but all non-gravitational experiments measure the 
"Jordan-Fierz" or "physical metric" . The "Jordan-Fierz metric" is related to the "Einstein metric" via the conformal 
transformation, 

9fiu = A 2 {kp)g lfliu . (2) 

Hereafter, we denote by a tilde quantities in the "physical frame" and by an asterisk those in the "Einstein frame" . 
From the variation of the action S we get the field equations in the Einstein frame 

G*ni/ — SttG^Tx^i; -\- T^^ L ij , (3) 
n„(p = -4wG*a(ip)T*, (4) 

where T^l is the energy-momentum of the massless scalar field, i.e., 

T^l = 2 ((p,p<p, v - ^9*^9^^^^ (5) 

and T?" is the energy-momentum tensor in the Einstein frame which is related to the physical energy-momentum 
tensor as follows, 

rr = = a 6 (<p)f^ . (6) 
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The scalar quantities T» and a(<p) are denned as T* = T*^ — T* h g*u V and a((p) = dhiA(tp)/d(p. It is apparent that 
a(tp) is the only field-dependent function which couples the scalar field with matter, for a(tp) = the theory reduces 
to general relativity. Finally, the law of energy-momentum conservation V V T" = is transformed into 

V w T t ; = a(<p)T m V mii <p, (7) 

and we set tpo as the cosmological value of the scalar field at infinity. In this paper, we adopt the same form of 
conformal factor A(ip) as in Damour and Esposito-Farese which is 

%) = > ! , (8) 

i.e., a(ip) = where (3 is a real number. In the case (3 = this scalar-tensor theory reduces to general relativity, 
while "spontaneous scalarization" occurs for /3 < -4.35 [13 ■ 

We will model the neutron stars as self-gravitating perfect fluid of cold degenerate matter in equilibrium. Then the 
metric describing an unperturbed, non-rotating, spherically symmetric neutron star can be written as 

dsl = g^ v dx»dx v = -e 2 *dt 2 + e 2A dr 2 + r 2 (d9 2 + sin 2 9d<j, 2 ) (9) 

where e~ 2A = 1 — 2fj,(r)/r while for the calculation of the mass function /x(r) and the "potential" function $(r) the 
reader should refer to Paper I. Finally, the stellar matter is assumed to be a perfect fluid 

V = (p + p) U^Uv + Pg„u. (10) 

where is the four-velocity of the fluid, p is the total energy density, and P is the pressure in the physical frame. 



III. BASIC PERTURBATION EQUATIONS 



In this section we present the equations describing perturbations of the spacetime, scalar field, and fluid in a spher- 
ically symmetric background. The equations we provide describe the non-radial oscillations of spherically symmetric 
neutron stars in scalar-tensor theories. We assume, in the physical frame, using the Regge- Wheeler gauge [23j, the 
following form of the perturbed metric tensor 



h =h^+h {+) 

ll fJ,V <<<u,U ~ l ilV 1 



(11) 



where h\u) denotes the axial (or odd parity) part of metric perturbations 
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and hffi denotes the polar (or even parity) part of metric perturbations 
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The functions hoj m , hi t i m , -ffo,z m , H\^ m , i?2,im> and Ki m describing the spacetime perturbations have only radial 
and temporal dependence while Y"; m — Yi m (9, 4>) is the spherical harmonic function. 

Following the previous definitions the perturbed metric tensor h*^ in the Einstein frame has the form: 
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where SA = A[3tp5ip is the perturbation of the conformal factor A and Sip is the perturbation of the scalar held, where 
Sip is a function of t and r only. The above definition of h*uv will be used to derive the perturbation equations: in 
other words, we will work out the perturbations in the Einstein frame and we will transform back to the physical 
frame whenever we need it. 

By defining the new set of perturbation functions Hq, Hi, H^, K, hg, and h\ as follows 

ASA) , H 1M = ^H Um , (16) 

ASA) , K lm = (K lrn - 2ASA) , (17) 

h>x,lm = ~^hij m . (18) 

The perturbed metric h %liv , in the Einstein frame, is simplified considerably and reduced to the "standard" Regge- 
Wheeler form of a perturbed spherical metric. We should notice that the scalar perturbations 8 A are linked only with 
the polar perturbation functions H$i rn , H\^ m , i?2,im, and K\ m . The axial perturbation functions hojm and h\^i m are 
only affected by the contribution of the scalar field to the background. 

The perturbation equations will be derived by taking the variation of Equations © and Q 

SG^ V = 8ttG*«V + ST^l, (19) 
U,Sp = -AttGJ [a((p)T m ] , (20) 

The various components of ST^l are expressed as linear combinations of Sip and h^ u . In the Einstein frame, the 
perturbed energy-momentum tensor describing the matter fields ST*^, is some linear combination of the velocity 
variation of the fluid 8U l ~ (WYi m ,VdeYi m — u sin -1 9d<pYi m , Vd<pYi m + u sin 9deYi m ) , the density and pressure vari- 
ations (Sp and SP) together with the variation of the scalar field (Sip or SA) and the metric perturbation h*^. The 
explicit form of the energy-momentum tensor is given in Appendix [XI 

The linearized Einstein equations (|19|l will be written as follows. From the tt, tr, rr components and the sum of 
the components 09 and (fxi>, we get 

EA ( ^™ = (/ = 0to3), (21) 

l.m 

where the four expressions Ar£ = are given in Appendix 151 They contain combinations of H a , Hi, H 2 , K, W , SP, 
Sp, Sp and their temporal and spatial derivatives. It is worth noticing that all four equations above are descibing only 
polar perturbations. In a similar way, from the tO, t<p, r6, and r<j) components, we get four more equations 



E {4^m + =0 (J — 0, 1), (22) 

l,m ^ ' 

E {&eY lm - 42^- e ^Yi m ) = (J = 0, 1). (23) 

l,m 

Here the expressions are some linear combinations of polar perturbation functions while on the other hand, the 

expression (3^ is a combination of only axial perturbation functions (see Appendix [BJi . 

Furthermore, from the Ocj) component and the subtraction of 90 and (fxj) components, we get two more equations 

J2{ s irnXi m -t lm sm9Wi m } = 0, (24) 
E {tlrnXim + &lm sin 6W im } = 0, (25) 

I .in 

where si m and ti m describe polar and axial type perturbations respectively while Xi m and Wi m are denned as 

X lm = 28^ (d e - Y lm and W lm = (d 2 g - ^-d - -^dl) Y lm . (26) 

^\ sinf / V sine* sur6» v 
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Taking the product of Equations (|21|l - l|25|) with Y/ m , integrating over the solid angle and looking at components 
with fixed values of I and m, we get ten partial differential equations in the variables t and r: 

4m=°> °t£=°> s lm =0, (1 = - 3 and J = 0,1) (27) 

/C )= 0, t/m=Q (-/ = 0,1). (28) 

Equations (|27|l describe the polar perturbations, and Equations (|28|l describe the axial perturbations. The analytic 
expressions for Equations i-e., Eqs. I|B7|I . (|B8|I and (jB10|) . do not couple to the perturbations of the scalar field 
5A or Sip. Therefore, the perturbed scalar field is coupled only to the polar perturbations. 

A. Axial perturbations 

It is quite easy to derive a wave equations for the axial perturbations by combining equations (|B8|) and (|B10(I 

X - e*~ A (e*- A X')' + e 2 * (^-^ - ^ + 4^G* (p - P) A^j X = 0, (29) 

where we introduce the new function X(t,r) defined as hi = e A ~*JTr. This equation does not couple to the scalar 
field perturbations, as we have mentioned earlier, and the effects of the scaler field will enter only via the background 
terms. Thus for (3 = 0, i.e., A = 1, it reduces to the standard wave equation describing axial perturbations [24| . 
Finally, from equation l|B8|l we get the following relation, 

u = -e- 2 */io, (30) 

which suggests that there are no axial oscillatory fluid motions, i.e., they have zero frequency and represent stationary 
currents. Thus axial perturbations are described only by a single wave equation l|29(l which does not couple neither 
to polar fluid or spacetime perturbations nor to the perturbed scalar field and it can be studied independently. 



B. Polar perturbations 



The equations describing polar perturbations can be simplified introducing a new set of perturbation functions. 
Introducing these functions we can reformulate the seven equations describing polar perturbations as a pair of wave 
equations and a constraint equation, using a procedure similar to Ref. 25]- The new metric perturbation functions 
will be 



e 2 * / - - \ 
F(t, r) = rK, and S(t, r) = — [H - Kj , 

while the fluid perturbations can be described by variation of the enthalpy function, i.e. 



(31) 
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The system of equations describing the polar perturbations can be reduced to the following pair of wave equations 
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From the tt component of the perturbed Einstein equations we get the Hamiltonian constraint: 
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Furthermore, for the special form of the conformal factor, i.e., A = exp(/3p 2 /2) or a — (3<p, from Equation (|2U[) we 
obtain a wave equation for the perturbed scalar field Sip: 



,2*-2Ar,„/ | 2* 



V + 



+ 4e -2A^2 _ 47rG>tA 4 ^ _ 3^ ( 4a 2 + ^ 



5yi 



r 2 e -2A^3 



2/^ 



inG^rA 4 {2r#P - a - 3P) | 



5 



2* 



2A lTl 3 , ^ +47rGU 4 | 2 ^p_ « ( -_ 3 p 



)} 



F - AnG^ae 2 ® (p + P 



- 3 ) H, (36) 



where G 2 = dP/dp. 

Finally, by linearizing equation 0, i.e., the energy- momentum conservation equations, we get a system of evolution 
equations for the components of the perturbed velocity and the density perturbation 
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Concluding, the dynamics of the polar perturbation is described by the system of evolution equations (|33|l . i|34|) 
and (|36l) together with the constraint equation (|35ll . The rest of the functions describing the spacetime and fluid 
perturbations will be computed by taking proper combinations of F, S, H and dp. 



IV. SPACETIME PERTURBATIONS IN SCALAR-TENSOR GRAVITY 



In Paper I we have studied stellar perturbations in scalar-tensor theories of gravity freezing the spacetime and scalar 
field perturbations. This is the so-called Cowling approximation, in which we only consider perturbations of the fluid. 
In practice we worked with a system of equations similar to l|3*7|) - (I51?!) . but setting Hi = H = H 2 = K = Sp = 0. 
Even under this approximation, spontaneous scalarization has a remarkable effect on the oscillation spectra of the / 
and p modes. Based on this observation we proposed in Paper I that a successful detection of gravitational waves 
from oscillating stars will provide us with a tool to constrain the phenomenon of 'spontaneous scalarization'. 

The quasinormal modes of the fluid perturbations described in Paper I will be affected by the coupling to the 
spacetime and scalar field perturbations. This is an interesting problem on its own. However, since we have already 
shown that the effect of spontaneous scalarization is quite strong when we limit consideration to perturbations of 
the fluid, in this paper we will examine the effect of the scalar field on the quasinormal modes describing the pure 
spacetime oscillations, i.e., the so-called w modes [2(j. The w modes are similar to quasinormal modes of black 
holes. They have higher frequencies and shorter damping times that the / modes, typical frequencies being of order 
7—12 kHz and damping times of order of 0.1 ms. These properties of the w modes are common to axial and polar 
perturbations |17| . In the case of polar perturbations the w modes are associated to small fluid motions, while in 
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the axial case there is no coupling with the fluid at all. This is the reason why here we choose to study the effect of 
the scalar field considering only the axial perturbations, described by the single wave equation (|2§|> . We expect the 
effect of the scalar field on the axial and polar w modes to be of the same order of magnitude. It should be pointed 
out here that, according to recent collapse calculations jS^, the w modes are significantly excited. This adds further 
motivation to our study of the effects of scalar fields on w mode oscillations. 

The equations needed to construct the equilibrium stellar configurations as well as the equations of state (EOS) 
used are described in Paper I. In Paper I we also considered cases where the asymptotic value of the scalar field 
(pa =/= 0, here, for simplicity, we only deal with scalar fields with ipo = 0. 

To compute the quasinormal frequencies of the axial w modes we will use two different techniques. In the first 
approach we carry out time evolutions of Equation l|29l) and Fourier transform the signal at infinity; in the second 
approach we assume a harmonic time dependence of the perturbations, and the corresponding boundary value problem. 



A. Evolving the axial perturbation equation 



The time evolution of the 1+1 equation (|29|) is rather simple to obtain. We set some arbitrary initial data (for 
example a Gaussian pulse) in equation (|29() and evolve these data for some time. Then we compute the oscillation 
frequencies by taking the Fourier transform of the signal emitted at infinity. 

In Figure ^ we show the waveforms observed at a distance of about 300km from a neutron star with Arnowit- 
Deser-Misner (ADM) mass 1.4M . It is noticable that the arrival time of wave for different values of (3 is not the 
same, because the effective potential due to central neutron star changes as a function of (3. In this figure, we can 
see that the waveforms for (3 = and for (3 = —4 are identical. This result can be understood as follows. For the 
axial perturbation, the gravitational wave is not coupled with the perturbation of the matter and the scalar field. So 
the presence of the scalar field is realized only due to the modified background. On the other hand, with ipo = 0, the 
central value of scalar field ip c is zero for any (3 > —4.35 (see 0]). Thus for values of (3 > —4.35 the effect of the 
scalar field is insignificant and it will affect the results only for (3 < —4.35. 
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FIG. 1: Waveforms describing the evolution of X for EOS A (left panel) and EOS II (right panel) generated for stellar models 
with Madm = 1.4M for f3 — 0, —4, —6 and —8. The cases j3 = and /3 = — 4 are indistinguishable. 



B. Boundary value method 



Our second approach to calculate the quasinormal frequencies of the axial w modes is more involved than simple 
time evolutions. However, using time evolutions we can only identify those modes that are significantly excited by 
a certain set of initial data. For example, using time evolutions it is not easy to identify quasinormal modes that 
damp out very fast. The approach we present here allows us to calculate both slowly and highly damped quasinormal 
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modes. We Fourier-expand the wave equation (|29|l as X(t, r) = X(r) exp(iujt) and get 
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(40) 



In this way we obtain an eigenvalue problem: the complex quasinormal modes uj can be obtained imposing appropriate 
boundary conditions. In our case, the boundary conditions are that X(r) should be regular at the stellar center and 
that there are no incoming waves at infinity. Inside the star we can just integrate the above differential equation; 
outside the star we use appropriate asymptotic expansions to ensure that there is no incoming radiation. Here 
we adopt a variant of Leaver's continued fraction method [28j. that has been originally used for the calculation of 
quasinormal modes of black holes. The procedure is described in detail in Appendix IU1 





Re(u>M ADM ) 



Re(ioM ADM ) 



FIG. 2: The frequencies of the axial w modes for 1 = 2, here we have considered stellar models with Madm = IAMq. The 
left panel corresponds to EOS A and the right panel to EOS II. We show both the "ordinary" w modes (wi, W2, W3, ...) and 
the wn modes (upper left corner in the plot). Open diamonds, squares and triangles correspond to the wu modes with (3 — 0, 
f3 = —6 and (3 = —8, respectively. Similarly, filled symbols refer to the fundamental w mode and its overtones. 



C. Results 



The results of the two methods we described agree very well, providing a good consistency check on our calculations. 
In Figures [21 El and^J we present the eigenvalues. Our results suggest that the presence of a spontaneous scalarization 
can be inferred from the w modes emitted by a newly born, oscillating neutron star. 

In Figure |21 we show the eigenfrequencies of the w modes for neutron star models with Madm = IAMq. The plot 
is reminiscent of earlier calculation of these modes (see e.g. Figure 3 in 1231 h The modes that might be relevant for 
gravitatonal wave detectors are the lowest w modes [2^. The wu modes[3(1| damp out roughly twice as fast as the 
w modes, but having lower frequencies they could also be relevant for detection by Earth-based interferometers. The 
higher-frequency w modes (u^, ^3, W4, ...)are difficult, if not imp ossible to detect. 

In the study of w modes as a tool for asteroseismology 0, 0, 0, it has been suggested that a proper 
normalization for Re(u>) is to multiply it with the radius R of the star and to scale it as a function of the compactness 
M/R. This phenomenological argument has been recently verified analytically by Tsui and Leung |3l| . Introducing 
/ = i?e(w)/27r, it is clear that Rf scales linearly as function of the compactness M/R. This applies both to Wu and 
u>i modes (and even to the higher overtones). The linear relations that can be derived from Figure El arc 

/ui -mode (kHz) = i [ax - Pi-jj^J and fun -mode (kHz) = i fan + Ai^J , (41) 

where the constants a±, (3±, otu and (3\\ are listed in Table [I] 

Another reason why it is harder to detect high-damped quasinormal modes such as the wu modes for compact 
stars is that the effective amplitude scales as the square root of the number of oscillations 22] . Typically we can 
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TABLE I: The coefficients for the fitting factors of equations (1411) . 



p 


oil Pi an 


Pn 





13.35 4.20 2.48 


3.32 


-6 


13.57 4.89 3.09 


1.88 


-8 


13.36 5.17 2.89 


1.50 



hardly observe more than 2 — 3 cycles for highly damped quasinormal modes of black-holes and for the Wn modes of 
compact stars. Spontaneous scalarization might help in this direction. Figure 0] shows that the damping time of the 
wu mode for stars with /3 < —4.35 is significantly longer than for typical stars in general relativity. The reason is 
that the presence of a scalar field increases the maximum mass of the stars and their compactness. Since the damping 
scales with compactness, the wu modes live considerably longer. On the contrary, the damping times of the w\ modes 
become shorter as the compactness increases (left panel in Figure 0J. 




-•-EOS A 
-■-EOS II 



0.12 0.15 0.18 0.21 



FIG. 3: The wi modes (left panel) and the wu modes (right panel) for EOS A and II and for values of /3 = 0, -4, -6 and -8. 



V. CONCLUSION 



We derived the equations describing stellar perturbations in scalar-tensor theories of gravity. The presence of a 
scalar field affects the equilibrium model, and consequently the oscillation spectrum. The scalar field perturbations 
couple with the polar perturbations of the spacetime and fluid, but they don't couple with the axial perturbations. 
Since the spacetime modes of polar and axial perturbations have the same qualitative behavior, we have chosen to 
study the effect of the scalar field on the axial perturbations. 

The results show that in the presence of spontaneous scalarization, a scalar field reduces the oscillation frequency 
of the wi modes by about 10% (i.e. by about 1kHz). The decrease in frequency for the wu modes is about 25% 
the frequency of (i.e., about 1.5 kHz). The effect on the damping time is even more pronounced. The damping of 
wu modes can decrease by as much as 30%, while it can increase by as much as 50% for the w\ modes. Detectors 
operating at these high frequencies are under development. Through a detection of the w mode spectrum, they could 
provide a unique proof for the existence of scalar fields with (3 < —4.35. 

A more detailed model of the effect of the scalar field on the oscillation spectra requires the inclusion of a larger 
set of equations of state. Another open problem is the study of polar oscillations, which couple directly to the scalar 
field. Work in these directions is in progress. 
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. (EOS A) 
= -6 (EOS A) 
= -8 (EOS A) 
= (EOS II) 
= -6 (EOS II) 
- 8 (EOS II) 



Wj mode 



S 

Q 0.42 
3 » .„ 



w n mode 



o a (EOS A) 

□ = -6 (EOS A) 

a !i = -8 (EOS A) 

• ,8 = (EOS II) 

■ .3 = -6 (EOS II) 

■ = -8 (EOS II) 



12 14 16 18 20 22 24 26 28 
M ADM /R 



12 14 16 18 20 22 24 26 28 

Madm/R 



FIG. 4: Dependence of the imaginary part of wi and ran mode on the stellar compactness Madm/R- The left panel is for the 
nil mode and the right for the ran mode. It is apparent that the imaginary part for wn mode can decrease by up to 30% in 
the presence of a scalar field. For the Wi modes, on the contrary, the damping time becomes shorter. 



Acknowledgments 

We acknowledge valuable comments by Emanuele Berti and Shijun Yoshida. This work was partially supported by 
a Grant for The 21st Century COE Program (Holistic Research and Education Center for Physics Self-Organization 
Systems) at Waseda University and the Pythagoras I research grant of GSRT. 



APPENDIX A: THE PERTURBED ENERGY MOMENTUM TENSOR 



In this Appendix we show the explicit form of the various components of the perturbed energy momentum tensor 
(of the fluid and of the scalar field) appearing in the perturbation equations. We will use primes for spatial derivatives 
and dots for temporal derivatives. For simplicity we will omit the subscript Im in the various perturbed quantities. 

The components of the perturbed energy momentum tensor T^l for the scalar field have the form 



01 *t4> 

ST iv) 



ST 



ST 



(v) 



ST 



„24>-2A 



2^>S(fi 



Y, 



Ira-) 



l 



ST ij) = e -2A^2 ha __ d4iY[ ^ 



(¥>) 

*ee 



sin 9 

-e- 2A V 2 h smed e Y lm , 
295<p?Y lm , 

2V6<pd g + e- 2A V 2 h 1 — 
sm( 

2*<5^-e" 2A * 2 /iisin! 

r_ 2 ^ + ^ 2 _ 

-2* V +(h 2 -K 



Yu 



r e 



r 2 e- 2A 



K 



). 



Im i 
2 



sin 2 m, 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 

(A6) 

(A7) 
(A8) 
(A9) 



In order to get the components of the perturbed energy momentum tensor for the fluid we define, in the physical 
frame, the variations of pressure 5P = 5PYi mi energy density Sp = SpY\ m and the components of the perturbed 
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4-velocity (in the physical frame) 



su f * = --* 



2A 3 

SU r = 4e 
A 



e HnYi 



*-2A 



WY lm , 



SU V 



Ar 2 



e* Vd e Y lri 



sin 9 



d&Yir. 



1 



Ar 2 sin 2 



-e* (Vd^Y lrn +usin9deYi m ) 



(A10) 
(All) 
(A12) 
(A13) 



where the perturbation functions <5P, Sp, W, V, and u defined in the previous relations depend only on t and r. Using 
the above definition the components of the perturbed energy-momentum tensor ST,^ V take the form 



ST^ tt 


= A 4 e 2 * 




= -A A e™ 


ST*a 


= -A 4 e 2 * 




= -AV* 




= A 4 e 2A 



p + P ) W + e~ 20 pHi 



u + e 2 ® pho 
-2* ~ 



p + p) u^y im - aV* ( /> - /' ) » r- < ■-">//„ 

4P 



1 



smtf 
sin 9d e Yi m , 



SP + PH 2 + —8 A 



A 



suit* 

5T W = A 4 p/j lS in^r im , 



A 4 r 2 



ST, 



A 4 r 2 



4P 

SP + PK + —SA 
A 



4P 

SP + PK + —SA 
A 



Y, 



! iii ' 



(A14) 
(A15) 
(A16) 
(A17) 

(A18) 

(A19) 
(A20) 

(A21) 
(A22) 



APPENDIX B: THE COMPONENTS OF THE LINEARIZED EINSTEIN EQUATIONS 



Here we provide the explicit form of the various expressions used in equations equations l|21|) - (|25|l for the description 
of the perturbed Einstein equations (jH?|) . We have chosen to use the same notation as Kojima to facilitate 
comparison. 



A 



(0) 



K" + e 2A [ -re~ 2A V 2 ' ^ 3 



r 



2r 2 



t 2X ( + \ ~ 8nG,A 4 p ) H 2 - 2*V - 87rG,A 4 e 2A (^ + IpcuV;). 



(Bl) 



A)2 = -K + e zly ( ^re-' 2A ^> 2 + ^ - - + AnG.rA^P I K + -H 2 + 



-,2A 



1 



-2A lT ,2 



3/X 1 



2r 2 



Hi 



- 2¥ty> + 87rG*A 4 e 2 * + P) W, 

( 2 ) - _ „-2*+2Aj> + e 2A /" I re -2A^,2 _ £. 1 



(B2) 



(Z -l)(i + 2) 



2r 2 



e 2A K 



Kl + V r 2Afr p 2A ( 1 



8nG*A 4 Pj H 2 - 2*<V - 8nG*A 4 e 2A (sP + 4Pa5<p\ , 



(B3) 
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A?2 =K" -HZ- 



Ira 



II _ -2#+2A 



K - II 2 ) + 2e- 2 *H[ - e 2A ( ^re^ 2 + ^ + 4ttG* (2P - p) taA H 



.,2 A 



^2 - 4 + - + IttG^P) H' 2 + e 2A ( 4^G* (P~p]rA i + ^—H± 



K' 



2e 



-2<t>+2A 



p — r 



4irG*rA i p H x 



fg + l) p 2A/V 
^2- e ^ 



„(°) 



- A ' 16^G*,4 4 P + ) /./,-+ -Mm/ - l<; ;r (V..J ; f iX (dF--\Fn<K-) . 

^ - e 2A (#2 + 



-2A 



„2A 



4ttG„ P-,5 rA 



14 ' "i) A 



8^G^ 4 e 2 * (p + p) 1/, 



(B4) 
(B5) 
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> + Pj rA 4 +(l- M j * 2 



{27TG, (; 

^3 {z(Z + l)r - Ap + 8ttG* (p + p) r 3 A 4 - 2r 3 e~ 2A * 2 } h - 87rG*yl 4 e 2 * (p + p) u, 



/? C1) 
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~2 


e- 2 *( 


'Af, - £1) 


Sim 


1 

~2 


(Po- 


-H 2 ) , 




=e 







-2* I 

-e h - 



2r 2 



i {2/i + 4^G* (P - pj r 3 A 4 } fc X) 



where the equations are simplified by virtue of the equations (13) - (17) in Paper I. 



APPENDIX C: NUMERICAL TECHNIQUES 



(B6) 

(B7) 
(B8) 
(B9) 
(BIO) 



In this Appendix we present two numerical techniques to determine quasinormal modes. The first is the direct 
evolution of the time dependent, axial perturbation equation 12911 . In the second approach we assume a harmonic 
decomposition for the perturbation function X of the form X(r,t) = X(r)e iuJt , and solve the equation (|40() as the 
eigenvalue problem. We shall consider the equations (|40|l in the interior and the exterior of the star; then we will find 
the eigenvalues (quasinormal modes) by matching the interior and exterior solutions. 

The largest numerical error in the interior solution occurs at stellar surface, where the pressure is zero. In order to 
avoid this difficulty, we integrate the perturbation equation (|40|l from both sides, i.e. from the stellar center r = 
and from the stellar surface r = R. Then we match the solutions at some intermediate point, e.g., r = R/2 (see for 
example [1^,|3^|). In order to deal with the boundary condition at infinity we adopt the continued fraction method, 
originally used for black hole perturbations by Leaver [2^. To use this method we must know the forms of the 
coefficient in the perturbation equation as functions of 1/r. Because of the presence of a scalar field , we do not know 
the exact forms of these coefficients. Therefore we just use the asymptotic forms of the coefficients, and derive a 
five-term recurrence relation. We believe that the QNMs obtained using these asymptotic forms are accurate enough, 
because the difference between the value of p at the stellar surface and at infinity is not so large. 



1. Interior region of the star 

The numerical integration of Equation (|40|l inside the star will be split (for numerical reasons) into two parts. First, 
we will integrate Equation (|40|l from the center towards P/2 and then we will integrate from the surface towards the 
same point. The matching of the two solutions will provide a unique solution valid throughout the star. 

Near the center it can be shown that X has a behavior of the form 



X = X c r l+l (1 + 0(r 2 )) , 



(CI) 
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where X c is some arbitrary constant. Using this boundary condition (jClf) and by integrating equation (|40|l from r = 
to the matching point r = R/2, one can obtain the values of X(r) and X'(r). For convenience we represent the two 
functions X and X' in the vector form Y = (X, X') and we will call Yo(r) the solution in the range < r < R/2. The 
next step will be to integrate equation (|40|) from the stellar surface towards R/2 with a set of boundary conditions at 
r = R such as (X(R), X'(R)) = (1, 0) and (X(R), X'(R)) = (0, 1). In this way we get two independent solutions Y x (r) 
and Y2(r) corresponding to each one of the previous boundary conditions. Thus the solution of the perturbation 
equation 14UI) is 

Y(r)=Y (r), for < r < R/2, (C2) 
Y(r) = oYi(r) + 6Y 2 (r), for R/2<r< R, (C3) 

where a and b are some constant, which will be determined from the junction condition at r = R/2: 

Y (R/2) = aY 1 (R/2) + bY 2 (R/2). (C4) 

The determination of the two constants specifies uniquely the solution in the interior of the star for a given value of 
the frequency cu and the constant X c . At the stellar surface the values of X(R), X'(R) are simply X(R) — a and 
X'(R) = b. 



2. Exterior region of the star 



The functions describing the stellar background simplify considerably outside the star. This leads to a corresponding 
simplification of the wave equation 1)41)0. In the exterior, the equations describing the background reduces to 



Vl = \re-^\ (C5) 

$' = \rtf + 4e 2A , (C6) 
2 r z 

V? = *, (C7) 

*' = -4(^-M)e 2A *- (C8) 



Therefore the asymptotic form of the above background quantities are 



H = M ADM + !±+o(\), (C9) 



r \r 



q> = _MADM +Q fl_ ^ (C1()) 



r 



^ = ^o + ^ + o(^), (Cll) 

where fii = —lo\/2 and uja = — ^ADM^s/^s' 

The perturbation equation (|40|l in view of the above relations, outside the star, get the form 



X" + ^X' 



2 -2* K l + 1) . 6/X 

Cu e j 1 — 3" 



X = 0. (C12) 



which is similar (in the absence of a scalar field, identical) to the Regge- Wheeler [23 equation describing the axial 
perturbations in the exterior of a spherically symmetric spacetime (either a black- hole or a neutron star). Using as 
boundary values for the integration the values of X and X' at the surface given by the two relations X(R) = a and 
X'(R) — b one can integrate equation (|C12(I together with (|C5J| _ (|C8I) from the stellar surface towards infinite. The 
numerical integration will obviously stop at some distance r = r a , where we will have to match the numerical solution 
with the appropriate asymptotic boundary conditions (in this case, the absence of incoming radiation) . 

In order to find the asymptotic form of the solution of equation l|C12|) we can assume a solution of the form 



*M=(4-l) X>„(l-^) , (C13) 
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where M = Madm ■ By substituting this form of the solution into the perturbation equation l|C12|) and taking the 
leading orders for \i and i.e., keeping only the terms up to order 1/r, from equations (|C9|) and HC10|) . we obtain a 
five-term recurrence relation for the expansion coefficients a n (n> 1), 

+ (3 n a n + InO-n-i + S n a n -2 + e n a n -3 = 0, (C14) 

where the coefficients of the recurrence relation are given by the following formulae 

Oi n — c n(n + 1), (C15) 

fin = don + c\n{n — 1), (C16) 

In = e + dx{n - 1) + c 2 (n - l)(n - 2), (C17) 

S n = ei + d 2 (n - 2) + c 3 (n - 2)(n - 3), (C18) 

e» = e 2 + d 3 (n - 3) + c 4 (n - 3)(n - 4). (C19) 

The coefficients Cj, <i and are functions of the background quantities and have the form 

i 2 ^ 2 ^ room 

c = 1 5-, (C20) 

Cl = _ 2+ ^ + % (C21) 

— l-^-Hf. (C22) 

c 3 = ™ + ^, (C23) 

C4 = (C24) 



do = -2iwr a - 2 + + ^ + JZL + JZL C25 

f Q r a r a A r a z 

12M 8*w/ii 24wM/xi 18/xi 

«1 = ^ 5 (U^DJ 

r a r a r a z r a l 

da = lf 3 M + 2^ 1+ 12 ^ /il + ^V (C27) 
r a y r a r a J 

d 3 = (3 + 4iwAf) , (C28) 

eo = — t(t + 1) + 2w ^1 H 1 2 1 J> (C29) 

'a 'a 'a 'a 

ei = — f-3M - 4c 2 M Ml + + - $^ , (C30) 

r n . \ r„ r a 1 



2u 

e 2 = 



i (3 - UloM} . (C31) 



The first four terms of the recurrence relation (|C14I) a_2, a_i, ao ; and ai are provided by the values of X and X' at 
r = r n , i.e., 



a-2 = a,-i =0, a = — — -, and a x = — — - 



X'(r a ) + ^^X{r a ) 
. r a -2M 



(C32) 



where 



-2iMu 
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The five term recurrence relations have in principle four possible solutions. A high order recurrence relation can 
generally be reduced to a three term recurrence relation, in which case convergence criteria for the solution can 
be easily applied, and we can identify the solution describing only outgoing radiation |2ftj |. To obtain a three-term 
recurrence relation we define new coefficients a ni /3 n , 7„, and S n as 

a n = On, P„ = Pn, In = In for n = 1, 2 (C34) 

and for n > 3 

n n O^n—l^n « Pn — l^n ? r "In — l^n irfn^\ 

a n = a n) p n = Pn ) In — In 1 , o n = n j . (C35) 

On-l On-1 «n-l 

The original five-term recurrence relation becomes a four-term recurrence relation: 

+ $ n a n + %a n -i + 6 n a n -2 = 0. (C36) 

Note that for the case of /xi = 0, that is the case of the standard neutron star obtained by Einstein's theory for gravity 
([3 = 0), the recurrence relation for a n has four terms pll l34l l35| . 

The final step will be to define another set of coefficients a n , f3 n , and 7 n : 

ai=di, j3 1 =0i, 7i =71, (C37) 

and for n > 2 

5 a Qn — lSn , ^ Pn—X^n / noo \ 

a„=o n , Pn = Pn ~ and j„ = j n : , (C38) 

7n-l 7n-l 

The four-term recurrence relation l|C36(l is thus reduced to a three-term relation of the form 

a n a n +i + (ia n + 7„a„_i = 0. (C39) 

Using this three-term recurrence relation, the boundary condition can be expressed as a continued fraction relation 
between a n , f3 n , and j n : 



a o Pi- P 2 ~ P: 
that can be rewritten as 



"■' ' " J <! " 2 "' : (C40) 



- a 7i «i72 a273 _ ^ fnAYS 
= p - = (C41) 

where /3o = ai/ao ; «o = — 1- The eigenfrequency w of a quasinormal mode can be obtained solving the equation 

/M = o. 
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